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Abstract 

To derive loop quantum cosmology from loop quantum gravity, I apply the model given in [48j to 
a system with coupled gravitational and matter fields. The matter sector consists of a scalar field 
(j) serving as a cosmological clock, and other fields {-0} providing physical spatial coordinates and 
frames. The physical Hilbert space of the model is constructed from the kinematical Hilbert space 
of loop quantum gravity, and the local observables in the physical Hilbert space are constructed 
using the matter coordinates and frames. A specific coherent physical state is then chosen, whose 
expectation values of the local observables give rise to homogeneous, isotropic and spatially flat 
gravitational and fields at a late clock time. The equations governing these fields may be derived 
using the symmetry of the physical Hilbert space. When the matter back reactions from are 
negligible, the result gives a specific loop quantum cosmological model in the 0{hP) approximation, 
with calculable higher order corrections. 

1 Introduction 

Loop quantum cosmology refers to a set of quantum cosmological models [Ml ESI [23 1211 HS] incor- 
porating key features of loop quantum gravity |35l EH |37| . As a candidate fundamental theory of 
gravity, loop quantum gravity has a rigorously defined kinematical Hilbert space describing quan- 
tum spatial geometry that is discretized in Planck scale. Building on this kinematical foundation, 
theorists are now tackling the challenge of obtaining the dynamics. By introducing the key features 
of loop quantum gravity into symmetric cosmological models, loop quantum cosmology aims to 
explore the effects of the quantum geometry in cosmic evolution. 

Loop quantum cosmology bypasses much of the complexity in loop quantum gravity because 
of symmetry reductions, and its dynamics is well-defined and extensively studied. However, the 
simplifications come with ambiguities: in the same physical setting, different symmetry reduction 
procedures can lead to distinct loop quantum cosmological models. Thus, several models in loop 
quantum cosmology have been proposed based on different physical considerations. However, for an 
isotropic, homogeneous and k > universe with massless scalar matter fields, several robust features 
prevail in the family of models [231 [261 123 (23 [22]: 1) they reproduce classical FRW cosmology 
in their large-scale semi-classical limits; 2) the semi-classical limits deviate from FRW cosmology 
near the initial singularity, replacing it with a well-behaved bounce of the scale factor [33l [33]; 
3) with a proper coherent state, an in-built slow-roll inflation phase occurs following the bounce 
|3H I32j. These results are obtained through both numerical and analytical approaches to different 
degrees of approximation. Moreover, the study of the cases with either negative k, inhomogeneity 
or anisotropy |24^ H5l |2T1 l22l [23] also finds a bounce in Bianchi I, H and IX anisotropic models. 

With its rich implications, it is essential for loop quantum cosmology to find its root in loop 
quantum gravity, so that the implications can be attributed to fundamental principles. On the 
other hand, people are searching for observable predictions in loop quantum gravity, hoping to 
derive loop quantum cosmology. 
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Aiming at these goals, I apply the model proposed in [48j to a cosmological system with A = 0. 
The matter sector of the system consists of a chargeless and massless scalar field (f) and other 
matter fields denoted as These matter fields live in the quantum geometry of space given 

by the gravitational sector of the system. The model describes this system using the kinematical 
Hilbert space of loop quantum gravity, called knot space. 

Based on knot space, the physical Hilbert space H of the model is built under concrete as- 
sumptions. Local observables in HI for gravitational and 4> fields are constructed, using spacetime 
coordinates provided by the matter fields [381 ES] • Specifically, the (p field serves as the cosmological 
clock, and the fields in {ip} provide the spatial coordinates and frames. For a suitable coherent 
state in H, the expectation values of these observables give rise to emergent classical gravitational 
and (p fields in the spacetime manifold. In our context, we will choose a coherent state that gives 
homogeneous, isotropic and spatially flat emergent fields at an early clock time. The cosmic evo- 
lution of the emergent fields will be calculated using the symmetry of H. When the matter back 
reactions from sector are negligible, the result conforms with a specific loop quantum cosmo- 
logical model at order 0{hP). To the approximation, it recovers FRW cosmology in large scales 
and resolves initial singularity by the big bounce. 

2 Kinematic Hilbert Space of Loop Quantum Gravity 

This section is a brief introduction to the kinematics of loop quantum gravity including gravita- 
tional and matter sectors. The kinematical Hilbert space, called knot space, gives a background 
independent description of the spatial quantum geometry and the resident matter fields. 

2.1 Ashtekar Variables for Gravitational and Matter Fields 

Loop quantum gravity is based on a canonical general relativity in the Ashtekar formalism \15\ [T6l 
[l7j . This formalism describes gravitational fields in a form similar to that of matter gauge fields. 

Traditionally, canonical general relativity uses spatial metric and extrinsic curvature defined in 
the spatial manifold M as phase space variables |15 1 [T6 l [T7]. As an alternative, triad fields consisting 
of three orthonormal vector fields {e^(x)}(x = (x, y, z); a = x,y, z;i = 1, 2, 3) can be used in place 
of the spatial metric. In the basis of e"(x) and its inverse ej,(x), the spatial Levi-Civita connection 
and extrinsic curvature take the forms rjj(x) and -R'*(x). Ashtekar variables (^Jj(x), £'?(x)) are 
related to {ef (x) , K'l^{x)) through a canonical transformation by |16)|17j: 

i^f(x) = det(e)e?(x); A^x) ^ rl{x) + jK{x) (2.1) 

where the real number 7 is called the Immirzi parameter. By construction, the fields E^(x) are 
densitized triad fields, and the fields ^^(x) are SO{3) gauge fields. The variables have the non- 
vanishing Poisson brackets: 

{Aj.(x),4(y)} = 87r{G/c')j5y^5{x,y) (2.2) 

where G is Newton's constant. Note that one may replace the 50(3) symmetry group with SU{2) 
in this formalism, since SU{2) and 50(3) share the same Lie algebra. 

Matter fields with a gauge group G can be included in this formalism [181 Ull ISO]- Using 
the triad basis, we describe the fermion, scalar and gauge fields by (^j(x), 7r|(x)), {(p\x) , Pi{x)) and 
AJj(x), E?(x)). Here i, i and i are respectively (gravitational spin ^) 5C/(2), Q and adjoint Q indices. 
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In this formalism, the action of general relativity possesses the following local symmetries: 1) 
time re-foliation invariance, which introduces the Hamiltonian constraint H{N) = 0; 2) spatial 
diffeomorphism invariance, which introduces the momentum constraint M{V) = 0; 3) local SU{2) 
invariance, which introduces the SU{2) Gauss constraint G{K) = 0; 4) local Q, invariance which 
introduces the Q Gauss constraint GiX) = 0. Each constraint is a functional of the dynamical fields 
and the Lagrangian multipliers -/V(x), F"(x), A*(x), A'(x) (the bars indicate their non-dynamical 
nature). The first three constraints result from spacetime diffeomorphism symmetry, and they 
consist of pure gravitational and matter terms: 

H{N) = Hg{N) + Hm{N)- G{k) = Gg{K) + G^{k)- M{V) = Mg{V) + M„,{V) (2.3) 

Here, we give the explicit forms of the pure gravitational terms, denoting k = 87r(G/c^): 

Hg{N) ^ d'xNix)Hgix) ^ (2ac)-i d'^Nix)^^ [e'^F^, + 2(1 + 7')i^f„^^] (x) 

Gg{K)= I d^xK\x)Gg,,{x) = {^K)-^ I d'xA^fdaEt + eij'AiEj^^ix) (2.4) 

JM JM ^ ' 

Mg{V)= f d'xV''{x)Mg,a{^) = {jK)-' [ d^xV'^{x)(EfFl,-K{l-^-')KiGgAix) 

JM JM ^ ' 

The explicit forms of the matter terms depend on the matter fields, which we would specify later. 

The three constraints (2.4) form a closed algebra with structure functionals independent of the 
exphcit forms of the matter terms. Denoting [A, A'] to be the SU{2) commutator, and setting 
[V, F']" = V'dbV"' - V'^dbV", the algebra is given by 

{G(A),G(A')} = 87r(G/c=^)7G([A,A']); {G(A),M(F)} = ^-^iG / d')^G{Cyk) 
{M{V),M{V')} = 87r(G/c^)7M([y,F']) 
{G(A), F(iV)} = 0; {M(F), H{N)} = %'K{G/c^)^H{CyN) 

2 /\VPaa AT Tpba ATI] 



{H{N),H{N')} = 8^{G/c')j {M{S) + G(5M„)) + ^^^7^ (^^^^|^^|p^ 



(2.5) 



where Cy denotes a Lie derivative and 



T^b rpai 

5« = {NdbN' - N'dbN) (2.6) 

2.2 Knot Space 

In QCD, the generalized electric flux and magnetic holonomy variables are powerful in capturing 
non-perturbative degrees of freedom of gluon fields. The Ashtekar formalism enables an analogous 
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treatment of gravitational fields. Gravitational holonomy and flux variables are signatures of loop 
quantum gravity [361 ESI EZ] , and they capture non-perturbative degrees of freedom of gravitational 
fields. The holonomy variable over an oriented path e C M (the bar indicates that e is embedded in 
the spatial manifold M) gives the parallel transport along the path by the connection fields Al(x). 
The flux variable over an oriented surface S C M gives the flux of -Ef (x) through the surface. 
Explicitly, we have: 

h^^){e)~^[A] = [VeKpJ^de'Ali^)T^^^f, F,{S) ^ f.E^dSa (2.7) 

where V denotes path ordering along e, and the S'C/(2)-valued gravitational holonomy is written 
in the spin (j) matrix representation. The matter variables compatible with the gravitational 
flux and holonomy variables are the following |T8l [HI |20] . The gauge fields AJj (x) are described 
by the Q holonomies h^'^(e)j in representations (i), and the E?(x) fields are described by the flux 

variables Fi(5). If v is a generic point in M, the ^y(x) fields are described by the irreducible tensors 
6^^\v) obtained from their Grassmann monomials of degree d, and the spinor momenta '7r|(x) are 
described by r]{v) = iO^'^^^^v)'^ . The (/>'(x) fields are described by h^^\v) = exp{(l^ {v)t^^^ ^) called 
point holonomies in representations (k), and the momenta Pi(x) are described by Pi{v). This new 
set of gravitational and matter variables are collectively called loop variables. The kinematical 
states of loop quantum gravity f36][35][37] are called knot states, and they are functionals of the 
configuration Ashtekar variables {j4^(x), AJj(x), ^*(x), (/)i(x)} via the corresponding loop variables 

{/i(j')(e),h«(e),0W(t;),/iW(?;)}. 

A knot state is given by an SU{2) x Q invariant product of loop variables defined on a graph 
in M, so it solves the Gauss constraints. Further, the knot state is only sensitive to the embedding 
of the graph up to a spatial diffeomorphism n G diffM, so it also solves the momentum constraint. 
Define an embedded graph 7 in M to consist of A^e smooth oriented paths {cj}, called edges, 
meeting at most at their end points {vn}, called nodes (the bars again indicate that 7 is embedded 
in M). Carrying loop variables, an embedded colored graph f is defined by: 1) an embedded graph 
7; 2) an SU{2) spin representation ji and a Q group representation jj assigned to each edge; 3) 
generalized SU{2) and Q Clebsch-Gordan coefficients (intertwiners) i„ and i„, a point holonomy 
representation k„ and a Grassmann monomial degree dn assigned to each node. The assignment of 
ji, jj, in, in! kn and d„ gives a SU{2) x Q scalar functional Sp: 

( Af„ Af„ Af„ Ne Ne ^ 

Sf[A, A, e, 0] = i (g) i„ (g) i„ (g) e^'^") {vn) (g) /i^"^") {vn) (g) h^'^~^ (ei) (g) h^^) (e^) \ [A, A, C, </>] 

( n n n n i i ) 



(2. 



where /nu{...} denotes the SU{2) ® Q invariant contraction. An element /i G diffM drags f to 
r' = /if and transforms Sf to Sff- = Sfi. Note that for every 7, there is a subgroup diffM, ^ of 
diffM that leaves the graph invariant, maintaining the set of all the edges and their orientations. 
There is also a subgroup TdiffM,-y of diffM,^ that acts trivially on 7. Thus we have the graph 



symmetry group Gm,'^ = dijfM,-y/ TdiffM,-y- The group averaging operator 
embedding information from Sf, and it is defined as: 



diffh 



serves to erase the 



Sf 



Sf 



Ng 



E 



A/, 7 



E 



= S 



[f] (2.9) 



4 



where Ng^^ is the number of elements in Gm^^, and [F] is a colored graph obtained from the 
embedded colored graph f by erasing its exact embeddingj^ The result S[f] is a knot state. Such 
a state is determined by a colored graph [T], and is diffM invariant since Sjf] • fi = s^^p^ = s^py 
The inner products between knot states are given by (generalized) Ashtekar-Lewandowski measure 
|36j |35j [37] [20j ■ with which the set of all knot states gives an orthonormal basis: 

(s[f,]|s[f]) = 5[f/]jf] (2.10) 

This basis spans knot space K, the SU (2) x Q and diffM invariant kinematical Hilbert space of 
loop quantum gravity. 

Canonical quantization of loop variables leads to the operators of the form 0{Q), where d C M 
may be v, e, S or R. However, 0(0) does not preserve K since Cl is not dijfM invariant. To 
respect diffM symmetry, we now replace O by a dynamical object Q that assigns (l{T) C M to 
each r, such that (^{T) transforms together with f under diffM transformations. Setting {r^ep} 
to contain exactly one representative from every [T], one can construct by specifying i).(Trep) 
for every T^gp G {T^gp}. More formally, each dynamical object Q is a map 17 : T — y r2(r) d A4^ 
satisfying Q{fiT) = fi' ij,Q{T) for any ^ G diffM and some /i' G TdiffM, ^j.'y- Then, we can define diffM 
invariant operator 0(0) as: 

• 0(J]) = Sp • d(n(f ))]P,,^^ (2.11) 

For example, one can construct diffM invariant flux and holonomy operators using a dynamical 
surface S : F — t- S{T) C M and a dynamical path e : F — t- e(r) C M: 

5[r] • mS) = 5r • F,{S{T))^d^ff,n m ■ ^^'^^4 = • hl^) {e{f)p ..g^ (2.12) 

The two operators are respectively differential and multiplicative operators acting on S[f], and their 
SU{2) X Q invariant products give gravitational operators in K. Since Ef{x) determines the spatial 
metric, the spatial area and volume operators are made up of the flux operators [361 ESI [37] . Fi{S) 
as a differential operator acts on the gravitational sector of Sf in a way that satisfies the Leibniz 
rule. Specifically, when S intersects with 7 only at its node vi, we have: 

Ny Ne Ny Ne 

F^{S) • (g)i„(g)/i(^'»)(e,) = i{S,eMei',vi)J0i,) • (g) i„ (g) /i(^'')(e,) 

Ny Ny 

Me) • (g)in(g)/i(^''He^) = ihKjY^ [<5e-g-, • /j0'=)(efc)7v(^'=) -<5_ _-i . Tl^'=h(^''\ek)] (g) in (g) /i^^'^e^) 

n i k n i^k 



(2.13) 

where i{S, Ci') is +1 or —1 when Cj/ is above or below S in an infinitesimal neighborhood of vi, and 
is if otherwise; i{ei',vi) is +1 or —1 when vi is the source or target of Cj/. 

The area operator in loop quantum gravity is derived from regularizing the classical expression 
in terms of the flux variables. The resulting area operator of a dynamical surface 5 is [8] 

Sir] -As^Sr- ^JF,iS{f))F^{S{f))P,,ff^ (2.14) 
1 Note that when f' = f/i, we have [f] = [f']. 
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Remarkably, knot states are eigenstates of As for any dynamical surface S. Set p to be a point 
where 7 intersects transversely with S{T), and let ji^ represent the spin carried by the edge Cj^ 
containing the point p. One can then easily check that 

S[f] • As = Mll-y V^(V+l)^[f] (2.15) 

p 

where the Planck length Ip is defined by /p = h{G/c^). It is clear that the spectrum of area operators 
in K is discretized, built from multiples of small quanta of order of Ip. 

The volume operator in loop quantum gravity also results from the regularization of the classical 
expression using the flux variables, but with much more complicated technicalities. Also, there are 
a variety of different regularization approaches that lead to different operator forms [9l [10] . Here 
we will follow the one introduced in [lOJ. The volume operator given by this approach acts on S[P] 
with a set of flux operators over dynamical surfaces intersecting 7 only at the nodes. With the help 
of (2.13), the action of the volume operator for a dynamical region R boils down to [lOj : 

S[r]VR = Sf^ J2 \/\(ivJ^d^ffM ; = ^ sgn{ei,ej,ek)eP''''Jp{ei)Jq{ej)Jr{ek) 

(2.16) 

By (2.16) the volume operator for a region R acts only on the nodes of a knot state contained in R. 
Further, each node-wise operation on the spin network state is given by the sum of the triplets of 
the action (2.13) on the edges meeting at the node. Therefore, the volume operator changes only 
the intertwiners of the nodes of [T]. Consequently, there is a volume and area eigenbasis for knot 
space, consisting of states given by a definite 7 and ji. 

Each member of the eigenbasis for the area and volume operators has the following interpreta- 
tion: a certain quantum of volume is assigned to each node of the state, and a certain quantum of 
area is assigned to each edge of the state. Through this basis, K provides a concrete description 
of quantum spatial geometry - the networks carrying Planck sized units of areas and volumes, 
associated with their edges and nodes. 

The remaining constraint to be imposed for a physical Hilbert space is the Hamiltonian con- 
straint. Adhering to the polymer-like structure of knot states, the standard Hamiltonian constraint 
operator H{N)^Q^ = Hg^^{N) + H^^'^{N) [Ml [SO] is quantized from a regularized discrete ex- 
pression approximating H{N). In the discrete expression, the curvature factors in H(N) are approx- 
imated by holonomies along a certain set of tiny loops. The quantization then leads to H{N)^^^ 
that contains holonomy operators based on these loops. Therefore, the action of H[N)^'^^ on 
a knot state involves a change in the graph topology that adds the set of tiny loops to 7. 

Moreover, since a non-constant Lagrangian constraint N is not diffM invariant, H^^^{N) does not 
preserve K in general. With constant N the action of H^^^{N) preserves but it is intricate 
and changes the topology of the graphs. As a result, constructing a physical Hilbert space annihi- 
lated by H^^'^{N) is a major challenge, which is currently tackled by both canonical approaches 
[351 EZl iQl 12] and the path integral formalisms [371 ill HH] . 

3 The Model 

In our model, we will modify H^^^{N) into a graph topology preserving operator H{Np) defined 
on K. The greatly simplified setting will allow us to apply the group averaging method to construct 
the physical Hilbert space IH of the model, based on certain concrete assumptions. A set of local 
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Dirac observables in H will be obtained using the clock constructed from the (p field and spatial 
matter coordinates and frames from the {'0} sector. 

3.1 Modified Hamiltonian Constraint Operator and Physical Hilbert 
Space 

For each embedded graph 7, we label each of its nodes with an integer n and each of its edges 
connected to the node n by an integer pair (n, i) (the range of i depends on n in general) For a 
given 7, we denote its node nhy vZ, and the oriented path starting from vZ and overlapping exactly 
with its edge (n, i) by ■ (fig.la). To each pair (e^ j) we assign a minimal oriented closed path 
e^ij that lies in 7, containing the outgoing path e^^ and incoming path (e^j)~^ (fig-la). Using 
these labels, we define a set of dynamical nodes {vm}, satisfying Vm(r) = w^^-^-) such that n{m) is 
one-to-one. Corresponding to a given {vm}, we define a set of dynamical paths {em,j} satisfying 

(r) 

— ^n(m) ^^'^-'^ that jijt) is one-to-one. The dynamical closed path 
Cm.ij' is then determined by the outgoing em,i and incoming e^- dynamical paths. We also define 
a set of dynamical spatial points p = {pk}, where k ranges from to infinity, such that vZ = Pfc(r) 
holds for exactly one k for every f and nj^ Notice that there are infinitely many distinct sets of 
dynamical nodes, paths and spatial points satisfying the above. This ambiguity comes from the 
arbitrariness of identifying the nodes and edges between different knot states, due to the absence 
of a reference background in K. 

The modification from H^^^{N) into H{Np)~ a crucial step for the model- contains two el- 
ements [?8]. First, N is replaced by a dijjM invariant lapse function Np^pf^), which is a function 
of a set of dynamical spatial points p = {pk}- Recall that the embedded spatial point Pfc(r) e M 
transforms under diffM just as T does. That means the classical counterpart of Np{pk) is the field 
A^(x) that transforms as a scalar field. Therefore, in proper semi-classical limits, the operator 
H{Np) is expected to approximate the diffM invariant Hg{N) instead of Hg{N). Second, the tiny 
loops (fig. lb) that define the holonomy operators in H^^'^{N) are replaced by the new set of closed 
paths (fig.la) that are contained in the graphs of knot states. By using the corresponding new set 
of holonomy operators, H{Np) preserves the graph-topology of knot states. Since H(Np) preserves 
both diffM symmetry and the graph topology, it is an operator in any subspace Kj- C K with 
one specific graph topology T. For simplicity, the model's kinematical Hilbert space is set to be 
^Ttorus C ^ with the graph topology of a lattice torus Ttorus-, which has nodes and six edges 
connected to each of the nodes. Restricted to Kji^^^^, we have 1 < n < A''^, and 1 < i < 6 for t>„ 
and en,i , and ^ j will be a square closed path overlapping with exactly four edges in 7 

The modification applied to the gravitational term of H^'^'^{N) results to the new gravitational 

^Note that each edge has two labels since it contains two nodes. 

^The spatial manifold M contains uncountably many spatial points, but we use only countably infinite 
set {pk} in correspondence to the discrete structure of K. Also, the one-to-one correspondence between 
and Pfc(f ) is obvious had we taken {pfc(f )} for any f to be the set of all points in M. Here we define {pfc(r)} 
to be a countably infinite subset of M, while maintaining this natural condition. 
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£(1,2). 





^(1,3,2) 



Figure 1: The left figure depicts an embedded graph 7, with one of its nodes labeled by n = 1 and the 
four edges connected to this node labeled by The figures in (a) demonstrate the corresponding 

definitions ofv2, ej^j and e^^j, while the figures in (6) demonstrate the corresponding definitions 
of el and el . that are used in H^^^(N). 



term Hg{Np), which acts on (s[f]| G Kj-^^^^^ as (set h{e) = h^^/'^\e)): 

4(1 + 72 



Pfc(r)=f„ 



(f)) X] sgn (em,i(r),emj(r),em,fc(r)) 



i,j,k=l 



i(e-™,(f))f [h{emA^))j, [Hf^^^{l),V^r)]] h-\emA^))^j [He^,{r))l [Hf^^^{l),V^r) 



xh-\e^^k{r)l /i(e„,fc(r)^-,y(r) 



where we have 



^^piPk\pi,{f)=v,^(r)) sgn{em,ii^),em,j(X),e„^,k{T)) 

Vm i,j,k = l 



h{em,i,j{T)) - h iem,i,j{T)))_ [h (em,fc(r)) 



Here, the total volume operator V"(j=.^ is defined as 



/i(em,fc(r))j_,V(f) 



(3.1) 



(3.2) 



^{f) = Y Y sgn (e„i,i(r),emj(r),em,fc(r)) e^^''Jp(e„i(r))Jg(emj(r))Jr(e„,fc(r)) 



48 

I'm i,j,k=l 



(3.3) 
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The modification applied to the matter term in H^^'-'{N) [20j results to the new matter term 
Hm{Np), whose explicit form depends on the matter content. Similar to the pure gravitational 
term, it is also constructed from the loop operators in Kji^^^^, and acts on S Kji^^^^ as: 



/iW(z;„,i(f))j,pi(t;„(f)),^5i(t;„,,(f)) 



•G) 



(3.4) 

where Vn,i denotes the end node of en,i- 

Finally, the Hamiltonian constraint operator for our model is given by the self-adjoint sum of 
the gravitational and matter terms: 

H{Np) ^ \{H'g[Np) + KXNp)) + \[H'g{Np) + KXNp))^ 

In order to compare the model with loop quantum cosmology, we will choose an identical 
physical setting. In this setting, gravitational fields are fully coupled to the chargeless, massless 
real scalar field (/>, while all other matter fields with small back-reaction are collectively named as 
m' . The classical constraints in this case are: 



H{N) = Hg{N) + H^{N) + F|^}(iV) 



M V det E - 



+{2k) 



-1 



1 



/ Ai Vdet E 



6'^EtE';{da4>)db(t> + {2KfP^ (x) + i^{^}(iV) 



G(A) = G,(A) + G^(A) + G|^}(A) 

= i^^r^ j^^ d\A^ (^daE'^ + e^.'^AiE^) (x) + + G|^|(A) 
= G,,4A) + G{^}(A) 



M{V) = MgiV) + M^V) + M|^,|(y) 

= (7k)-^ / d^xF"(x) (eIf!^ - k(1 - ^')KGg, ) (x) 



M 



M 



(3.5) 



Note that the functionals Hg^^{N)^ Gg^^[h.) and Mg^^{V) are the constraints for the subsystem of 
only gravitational and fields, and themselves satisfy the algebra (2.5). 
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The loop operators for (p field are the point holonomy opreators h^'^\v) = exp(iA0(?;)) with 
arbitrary real A, and the conjugate momenta operators p{v). Their non- vanishing commutators are 
given by 

[U^\v),p{v')] = ifi5^^^>{i\)U^\v) (3.6) 

Corresponding to (3.5), we have H'^{Np) = H'^{Np) + H'^^^{Np). The modified operator H'^{Np) 
formed by quantizing H(f){N) is given by [20]: 

_Ql2 1 

S^^^%{Np) ^ S. _ _ . (^^^7)-V E^p(P^lp-.(r)=.^(f))^^(^>n(r))^.(t^rn(r))IPd.// 



Pfe(r)=i)„{r))p(^^m(r))p(vm(r))<j ^sgn (e„,i(r), emj(r), e^,fc(r)) 



(/l(e„,i(f)))^,y5(7;„(f)) 



/i-i(e„,,(f)) 



^ 1 

/i(e„^,fc(^)))_ ,y5(^;^(r)) 



Mem,,(r))) F5(^„(r)) 



(3.7) 



where 



^^(!;„(f)) = J^sgn (e™,,(f),e„j(f),e^,fc(f))A(^o)(r;„(f))-i [/i(^«)(e^,i(r)(i)) - A(^o)(r;„(f )) 



x(r«)^ (/i-He-^,,(r)) 



Mem,fc(f)))^F3(i;,„(f)) 



(3.8) 



with Ao a small real number whose exact value is unimportant in our context. The modified 
Hamiltonian constraint operator H{Np) in our setting is then given by: 



H{Np) = -iH',{Np) + H'M)) + ^{H'giNp) + H'^{Np))^ = Hg^Np) + -{H[^^ + H'l^^){Np) 

(3.9) 

Recall that our kinematical Hilbert space -f^Ttor^s is already SU{2) x Q and diffM invariant. To 
obtain the physical Hilbert space of the model, we still need to impose the remaining symmetry 
generated by {exp(ii?(A'p))} with arbitrary Np based on arbitrary;). These unitary operators form 
a faithful representation of a group G, that is 

oo 

n exp(ii?(iVpJ) = [U{g)] (3.10) 

k=\ ) ^Kr, ^ 

' torus 

where Np^ is an arbitrary lapse function based on an arbitrary pj^, and ~ K-ji^^^^ means that we 
identify two expressions if they give the same operator in Kjiorua ■ The physical Hilbert space of the 
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model is constructed using group averaging procedure under the assumptions: (1) the existence of 
the left and right invariant measure dg for G; (2) the operator P defined by 

¥= f dgU{g) (3.11) 
Jg 

maps (^1 G K-ji^^^^ into P]?/)) G X^^^^^. The two conditions hold in minisuperspace models [HI [1112]; 
but remain to be proven for the model. Under these assumptions, the inner product between any 
two states = Pj-i/'i) and |^2) = IPIV'2) may be defined as 

{^l\^2) = = (V'l|lP|V'2) (3.12) 

The physical Hilbert space of the model HI C -^-f^^,.^^ is the space spanned by {P|sjp])}. By 
construction, H is invariant under the action of {exp{iH (Np))} with arbitrary Np, and satisfies 
the modified Hamiltonian constraint. Each element in HI is a solution to the quantized Einstein 
equations, and therefore is a quantum state of spacetime with the matter fields. 



3.2 Local Dirac Observables 

The model obtains its local Dirac observables using clocks, spatial coordinates and frames given by 
the matter fields. For an arbitrary set of dynamical nodes {vm}, one may construct a set of self- 
adjoint, commuting matter operators consisting of scalar operators {(p^{vm), (p^ivm), (p'^ivm), (p^ivm)} 
diagonalized by the knot state basis of Kji^^.^,^, current operators {Vj{vm), Ujivm)} and conjugate 
current operators {V/ (vm) , Uj (vm)} {I = 1,2,3 for the vector currents; 7 = 1,2 for the spinor 
currents). The operator triplet {(j)^ (vm) , (p'^ (vm) , (p^ {vm)) = ^{vm) will serve as spatial coordinate 
operators, {Vj (vm) , Uj{vm)} and {V/ (vm) , Uj (vm)} will be spatial frame operators, and 4P{vm) 
will be the clock operator. In our setting, the cosmological clock field (t)^{vm) is chosen to be: 

^\v^) = (2ao)-M^^^°^(^^m) - h^-^'\vm)] (3.13) 

provided by the (f) field. 

Classical gravitational fields in the Ashtekar formalism are SU{2) tensors. In the model, the 
SU{2) invariant components of the fields are described relative to the matter spatial frames. Ex- 
plicitly, for any S[f] G i^Ttor^s we define 

S[r] • h(en,kYj = • ^(t;„,fc(r))f/i(e„,fc(r));;;7(t;„(f ))^}Pd,j/ 

In canonical general relativity with the matter spatial coordinate field <I?(x), we can obtain a diffM 
invariant and spatially local variable 0{X) by integrating det((?a*J'(x))(5(<I>(x) — X)O(x) over M. 
Analogously, the model uses a normalized Gaussian distribution 5'' with finite width e and defines: 

S[r]d{X) = 5r^det(A$(T;„(r)))5^(d(t;„(f)) -X)0(r;„(r))Prfi;^ 

n 

S[r]d'{exAx) = Sf^Yl det(Al>(?;„(f ))) det(Al>(t;„,i(f )))5^($(i;„(f )) - X)6^MvnA^)) - X - AX) 

n,i 

(3.15) 
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for any sjf ] G -f^Tiorus ' where the coordinate volume element operators are given by: 

It is crucial that the spatially local operators obtained in (3.7) do not depend on the choices of the 
dummy variables {vm} and {en,i}- The diffM invariant operators are localized by referring to the 
spatial matter coordinates. 

To specify time using the clock field (jP, we construct an operator IIt that maps a spacetime 
state I*) G H to the spatial slice state 11^1*) G ^Ttorus which = T. With cb{v) denoting an 
operator acting on v, the operator fix is defined by: 

r_ . . 1 . (3-16) 

where the symmetrization sym{...} in the ordering of n is applied. The resulting local Dirac 
observables in H are given by 

d{X, T) = PO(X)nr ; 6' {ex ax, T) = PO'(ex,Ax)nT (3.17) 

In the following, uo in (3.16) will serve as an index running over the matter coordinates and frames 
{(fp, <I>, f-Vi, f-V^, g-Uj, g-U^} contracted with the non-zero S'C/(2) test functions {/*, /j, 3*, ^j}. 
The variables thus contain information about the momenta of the reference matter fields. 

Finally, applying (3.14), (3.15) and (3.17) to the fields of concern wc obtain the relevant localized 
observables. In the following, wc will give the dynamics of Jiex.AX ,T)j, h{ex,AX ,T)j, p{X,T), 

h^^\X, T), in a certain "physical gauge" of matter coordinates and frames that is fixed by C'uiX, T). 



3.3 Conditions on Matter Coordinates and Frames 

Our next step is to choose a proper physical state |^') G BI to derive semi-classical limits from the 
model. In order to give sensible descriptions, the localized observables must refer to well-behaved 
matter coordinates and frames. So there are conditions to be imposed on the matter sector of |^'), 
which provides the matter coordinates and frames. 

For the clock, we require any two spatial slices of |^) with (fp = Ti and cj)^ = T2 to be related 
by causal dynamics, so that either one of them can be used to reconstruct the spacetime. In the 
model, this condition is imposed as: 

FUt.I'^') + 0{h) =FflT2\'^) + Oih) = 1$) (3.18) 

When (3.18) is satisfied, we have: 

Fdi{X)d2{X)UTM 

= Fdi{x)d2{x)iiT,Ffinm + o{h) 
= Fdi{x)iiT,Fd2{x)iiT,\^) + o{h) 
= di{x,Ti)d2{x,Ti)\^) + o{h) 
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More generally, denoting by (p'{Oi, O2, Oat) a function of linear operators in a Hilbert space, 
we have 



F^'{d,(x),d2{x),....,dN{x))tiTM 
= ^'{di{x,Ti),d2{x,Ti),....,dN{x,Ti))\^) + o{h) 

On each spatial slice IItjI^), we will also impose spatial coordinate conditions to properly 
describe gravitational and (j) fields. Denote an arbitrary function of the gravitational and (f> field 
operators based on {e^^j} and {fn} as: 

(p = ¥'(J(en,j)j,/i(e„,i)i,/i^(e„,i)t,p(i;„),^W(i;„),/iWt(^;„)) 

We require that there exist {v*} and {e* j} such that: 

= ¥'(J(ex„,Ax„,j7, h{ex^,AXr.,S hHex^AXr^Spi^n), h(^\x^), hW\Xn))fln 1^) + 0{h) 

(3.21) 

with a certain set of values {Xn} and {AXm,j} satisfying {X^ + AX^j} = {Xn}- Here, {f*} 
represents the physical nodes at clock time Ti for an observer using the spatial matter coordinates 
$ = {(p^ , (j)'^ , (p^) . Naturally, there is also a set of physical spatial points p* = {p^} agreeing with 
{f *}, such that Prni^) = v^i(^) for m < Ny. Condition (3.21) says that each physical node acquires 
a matter spatial coordinate vahic, and that the matter frames are physically orthonormal to each 
other at each physical node. The spatial coordinate condition can now be imposed on the map 
(f*,e* j) — >■ (Xn, A.Xn,i)- Analogous to the coordinate maps on a torus manifold, the map should 
appear smooth in large scales at most of the physical nodes except at some {v^^} C {w*} where 
coordinate singularities occur. We demand |AX„^j| < d for {Vn^}, where d is small enough 
that the map appears continuous at large scales. Also, for any m the set {AXm,i < d,} should define 
a parallelepiped in up to an error of 0{d) so the map appears smooth at large scales. Notice 
that once (3.11) is satisfied, the coordinate conditions can be achieved easily through a redefinition 
of coordinates 4 $'(!>). 

Using a lapse function satisfying N^{p^) = J\f{Xm) with an arbitrary function J\f{X), we 
can apply (3.21) to Hg^^^)- 

Hg,^{NfiUT, I*) ^ ^(A^(p;), J(e*,J,, /.(e;,)^, hHe:^,)],p{v:), U'\v*J, /^Wt«))n^ |*) 
= Hp,^(AA)nrJ*) + 0(/i) 



piXn),h^^\Xn),h^^^\Xn)) 

(3.22) 
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Copying the form of H'g{Ng,) shown in (3.1), n'g{N) is given by: 
VlgiN) ^ n'fiM) - 2(1 + 7^) g^, ^ Y.^{X) (^^.Ax, exAY, ^XAz) 



X 

J 



{h-\exAx))j [(/.(ex,Ax))_ , [Wf (1),V]J {h-\exAy)) j [(/i(ex,Ay)) _ , [wf (1),V 
(/i-'(ex,Az))^[(Mex,Az))' V 



(3.23) 



where HfiM) is 



U'^{M) = ^Q^2r^^j^f^^ Yl-^^-^^ H sgn(ex,AX,ex,Ay,ex,A^) [hiexAXAv) " /i~^(ex,AX,Ay)j 



h{exAz) ) - . V 



and 



V = XI X] (ex,AX, ex,Ay , ex.Az) e^'^^J (ex,Ax)pi (ex,Ay)Qi (ex,Az)i? 



(3.24) 



(3.25) 



Copying the form of H'^{N^) shown in (3.7), n'^{M) is given by: 



where 



-8V 



642 . . 2k 



A2 



XAA(X)^Q(X)^q(X) 



642 . 2^ 



(i;iK7)-'^XAr(X)]5(X)p(X) X sgn(ex 



272 • 482 

X (^"^(ex,Ay))^ 

X {h-\exAz))l {kexAz))''- W 



h{exAY))^,VHX) 



AX,Ay,AZ 

1 



^ (ex.Ay) , 
2 



Ay) _ ,vn^) 

M 



J2 sgn (ex,AX, ex,Ay, ex,Az) h^^\x)-' [h^^\x + AX) - h^^\x) 

AXAYAZ 

M ' 



:(r«)| (h-HexAv))'^ [(M(ex,AY))^ , W 
{h-\exAz)Y- \{HexAz)Y^,vi{X) 



(3.26) 



(3.27) 
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Further, the well-behaved matter coordinates and frames also lead to simple algebraic relations 
between the spatially local operators acting on (^IItiI*) = \ip). Prom now on we set a variable 
{X,AX) to take values in AX„^j)}. Also, define Q;(ey^Ay) as the operator localized from 

ot{^n,i) which is diagonalized by the basis {(s[r]| G ^Ttorusl- The state (sjpjl has eigenvalue +1 
if the embedded edge of F overlapping with e„^i(r) has the same orientation as 6^,1 (F) and has 
eigenvalue —1 if otherwise. The algebraic relations between the spatially local operators are the 
following: 

j\exAx)iW) = J{ex,Ax)i\<p) +0{h) 



J{ex,Ax)i,h\eY,AYYj Iv?) = (^x,y(^AX,Ay^^p7('r/) j ^Hey,Ay)^|</') 

-Sx,Y+AYS-AX,AYilll{T^)Rh\eY,AY)j\<p) + O(ip^) 



Hex,Ax)hkeY,AYYj\ \^) =0 + O{llh) 



hiexAx)h h\eY,AYYj\ \^) = + 0{llh) 



J {ex, ax) I, J (ey, ay) J \^) = 5x,y5AX,Ay«^p7e/J'^-^(ey,Ay)xa(ey,Ay)|<y?) + O^p^) (3.28) 



p{X)\cp)=p\X)\^) + 0{h) 



'p{X), /iW(y)l \^) = dx,YihiiX)h^^\YM + 0{h^ 



p(x),AWt(y)l 1^) = -5x,Yih{i\)h^'''^\Y)\^) + o{h^) 



h^^\x), h^^'\Y)\ |(^) = + Oih^) 



h^^\x), /i(^')t(y)] \ip) = + 0{ti 



\p{X),p{Y)]\^)=Q + 0{h^) 

This conditional algebra enables further calculations after the approximations (3.21) are made. 
Lastly, we also require 1^*) to give the momenta of the matter coordinates and frames: 



i>^{X,T{)\^) = u^{X,Ti)\^) + 0{h) 



(3.29) 
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The matter coordinates and frames specify a gauge for the gravitational sector. Locahy at the 
moment Ti, this gauge is characterized by (3.21) and (3.29). For this paper, we require the simple 
gauge conditionj^ 

u^oiX, T) / ; iy^{X, T) = Uf.Vj{X, T) = vj.yi{X, T) = Ug.Uj{X, T) = u-,pi{X, T) = (3.30) 
which gives a comoving frame in the cosmological setting that will follow later. 

3.4 Coherent States and Emergent Fields 

We have required the state |^') to satisfy the quantum coordinate conditions (3.21) and (3.29), such 
that the set of local observables { J(ex,AX, T)/, /i(ex,AX5 7") j} and {U--^\X, T),p{X, T)} would give 
meaningful descriptions around the moment Ti. Since our goal is to obtain semi-classical limits, 
we impose coherence conditions on the gravitational and (p field sector of |^') as: 

J{exAX,Ti)i\^) = [{W{exAX,Ti)m) |^) + 0(/2) 

(3.31) 

/iW(x,ri)|^) = ((^|/iW(x,ri)|^')) \^) + o{h) 
p{x,Ti)\^) = {{^\p{x,Ti)\^)) m + o{h) 

Since the clock cj)^ is built from cp field, the conditions imply: 

z>^o(X,ri)|^) = {{^\j)^o{X,Tim) m+0{h) (3.32) 

which specifies the value in (3.30). Additionally, we also want the expectation values to appear 
continuous in terms of the spatial coordinates for a semi-classical state. Therefore, for any two 
e* j(r) and ej^ ^(r) that share a common node and form a smooth path, we will impose (recall that 
\AXn,^\ < d for < {<J): 

Wiex„AX.,Ti)j\^) = {^\J{ex,^AX,,Ti)i\^) + 0{d) 
{^\hiexM,Tiyjm = mh{ex^AX„T,Yj\^) + 0{d) 

(3.33) 

mU^\Xn,Ti)\^) = {^\h^^\x^,Ti)\^) + 0{d) 
{Mp{Xn,Tim = {^\p{X^,T,m+0{d) 



^Note that the vaules of T) can be adjusted by clock-time dependent redefinitions of the matter 

coordinates and frames 
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Because of the algebraic relation (3.28), we expect the gravitational and (j) field sector's solutions 
to (3.32) and (3.33) to exist. These conditions say that |^') has sharply defined, approximately 
continuous values for the local observables at clock time Ti. Therefore, l^*) is expected to be 
semi-classical around that moment. 

To make contact with classical general relativity, the model maps the expectation values in 
(3.31) to the classical field values, using the matter coordinates as a common reference. For an 
explicit example, we will first pick a simple spatial matter coordinate system. Recall that the lattice 
torus Ttorus can be constructed by identifying the opposite boundary faces of a lattice rectangular 
prism l| C M^, which consists of the vertices {Vn = X^} and links {/j} (the bars indicate the 
embedding in M'^). Such construction naturally gives a coordinate map f * — )• t^i = Xn- Also, the 
approximated spatial coordinate space is naturally a rectangle region 7^ D 7| inside of . 

Once matter spatial coordinates are chosen, our model identifies every e* ^ disjoint from 
with an embedded path in 7^, under the guidance of the matter coordinate values. In our case, e* ^ 
is identified with the oriented path ex„,AX„ ^ that goes from the vertex Xn to the vertex X„ + AX„^j, 
which overlaps exactly with the link lii connecting the two vertices. Subsequently, we choose a cell 
decomposition dual to 7|, dividing 7^ into a set of cells {cx„} that are parallelepipeds up to errors 
of 0{d). Each cell cx„ uniquely contains a vertex Xn-, and the boundaries of {cx„} consist of a set 
of faces {si} whose each element Sj' intersects transversely with a unique link /j/. Suppose U' links 
Xn and Xn + ^Xn^ii we denote Sx^AXni ^ the oriented surface overlapping with Sj' and having 
the same orientation as ex„,AX„ ^ Then, model uses a fitting algorithm that maps the expectation 
values {{J{ex,/ix-,T)i) , {h{ex,AX7T)j)} and {{h^^\X,T)), {p{X,T))} to the values of the smooth 
fields {Ef{X, T),A( {X, T)] and {0(X, T),P{X, T)} defined in 7^. The fitting algorithm is required 
to obey the following rulesj^ 

/ Ef{T)dsa = {J{exAX,T)i) 
Pexp[/ Ai{T){Tj)de''\f = {h{exAX,T))f 

•'ex, AX 

(3.34) 

/ P{X\T)dX' ^{p{X,T)) 

J cx 

exp{iXcPiX,T))^{h^^\x,T)) 

In the rest of the paper, we will use the spatial coordinates with AXm^i £ {{^d, 0, 0), (0, ±d, 0), (0, 0, id)} 
for {^^nfc}' which makes cx„ a right cubical cell, and Sx,AX a square oriented surface. The 
choice of 7| C M'^, the cell decomposition {sj}, and the fitting algorithm described above are re- 
stricted but non-unique. However, any choice satisfying the restrictions gives a valid correspondence 
between \^) and the emergent fields. 

^ Note that such an algorithm is guaranteed to exist, since we are fitting the smooth fields with infinite 
degrees of freedom to the finitely many data points given by the expectation values of the local observables. 
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4 Physics in Homogeneous, Isotropic and Spatially Flat 
Sector 

The paper |4.8j has shown that the fields Ef{X,T) and A^{X,T) determined by l^') reproduce 
vacuum general relativity up to quantum gravitational corrections, when matter back reactions can 
be ignored. In the following a similar calculation will be done, assuming l^') gives homogeneous, 
isotropic and spatially flat emergent fields {Ef{X,Ti),A'^{X,Ti)} and {4>{X,Ti), P{X,Ti)} at a 
late initial time Ti. In comparison with [IH] the calculation will be done with two extensions: 
1) while the matter back reactions from {ip} will still be ignored, the full interaction between 
gravitational and (j) fields will be considered; 2) the quantum gravitational corrections will be 
evaluated to 0{hP). The calculation will give the 0{h^) effective equations governing the fields 
{Ef{X,T),A({X,T)} and {4>{X,T), P{X,T)}, which will be compared with the 0(^°) effective 
equations in different models of loop quantum cosmology. 



4.1 Emergent Constraints and Diffeomorphism Algebra 

For notational simplicity, we will denote the collections of the spatially local gravitational, (p held 
and a operators as {J,h,p,h^^\a}. One can evaluate the n-fold commutator by first applying 
(3.28): 



n-l 



(4.1) 



(^|]p$„( J, h, h^^\p, d,M)iiTi 1^) + o{h) 



where is given by carrying out all n commutators. Then we use (3.19), (3.31) and (3.34) to 
obtain: 



UgANn),.... [-Hg.^AAs), \ngAN2).ng,4,{Nl) 

{^\^n{J{Ti), h{Ti), AW(T^i),P(Ti), a{Ti),Mim + 0{h) 

<i>„((j(ri)), {h{Ti)), {h^^\Ti)), iPiTi)), (a(ri)),M) + oih) (4.2) 

{Hg^Nn), .... {HgAN3), {HgAN2),HgANi)}} ...} \ (T^) ^A^ {T^) ,4>{T,_) ,P {T,) ,N i=Ni 

+0{h) + 0{d^) 



Recall that Hg^Mi) is the Hamiltonian constraint for the subsystem of only gravitational and (j) 
fields. The {■0} sector serves as a background for this subsystem. Since the spatial coordinates are 
given by the {■0} sector, the lapse functions Mi{X) are Lagrangian multipliers for the subsystem. 
According to (2.5), the Poisson brackets, to all orders in n, with arbitrary A/i reproduce the full (off- 
shell) diffeomorphism algebra between Hg^N), Gg^K) and Mg^(^(y) in the semi-classical limit of 
\^), up to corrections of 0{h) + 0{d'^). 

Moreover, the symmetry H{Np)\^) = implies that: 



n-l 



H{N- 



■Afn 



F(<^),i?(<- 



HtiI^') = 



(4.3) 
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Separating the contributions involving H'^^^ or H'^^^ and denoting them as e{^}, we have: 



= 



n-l 



n—1 



n—1 



H{N^."),.... H{Np), H{Np),H{N; 



pV2 



nnl^-) 



nTi|*) + ew 
JItA^) + 0{h) + e{^y 



$„((J(ri)), Chin)), (/i(^)(Ti)), (|}(ri)), (a(ri)),AA,) + Oih) + e|^| 

{ifg,<^(iV„), .... {Hg^^{N3), {Hg^^{N2),Hg^4,{Nl)}} ...} Ua(y^)_^j/(r^)_^(T,^)_p(r^)^^__^. 

+0{h) + O(d^) + e 



(4.4) 



where we replace Hg_^{N'^''') with i-Lg,4,{Mn) using the spatial coordinate condition (3.21). The 
term e^^j. represents the matter back reaction from the {■)/'} sector, and we will use e^^j to denote 
generic matter back reactions from the {V'} sector in the following. According to (4.4), the emergent 
gravitational and (j) fields satisfy: 



'^9,'t>mE^^{T,),Ai{T,U{T,),P{T,) = ^ + 0{K) + 0{(f) + e{^} (4.5) 
MgA^)\E1{T,),Ai{T^U{Ti),P{T,) = ^ + ^ifi) + O(d^) + e{^} 



Therefore, the emergent gravitational and (j) fields as a subsystem are on-shell up to corrections 
0{h) + 0{d^) + e^^y 

4.2 Homogeneous, Isotropic and Spatially Flat Dynamics 

The coherent state \^) gives the dynamics of {Ef{X,T),Al{X,T),(t){X,T),P{X,T)}, through 
the corresponding {{J{exAX,T)i) , {h{ex,AX,T)^j), {h'^^\X,T)), {p{X,T))} at various T. We now 
evaluate this clock time dynamics in a symmetrical setting. 

Before the calculation, we need to set up a homogeneous, isotropic, and spatially flat ini- 
tial condition at some late clock time Ti. Since the (p^ field serves as the clock, it is clear that 
{(jP(X,Ti)) = Ti, and (3.13) implies that (p{X, Ti) is spatially constant. We require the coherent 
state 1^') to give homogeneous Ef{X,Ti), A^{X,Ti) and P{X,Ti), such that: 

Ef{X, Ti) = Ei5f ; Ai{X, Ti) = Aid^ ; P{X, Ti) = Pi ; (/.(X, Ti) = (/.i (4.6) 

where (5" identifies {x, y, z) with (1, 2, 3). Note that (4.6) satisfy (4.5) as required, and it also implies 
that 

v^{X,T^) = ui (4.7) 
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The dynamics of the emergent fields results from the symmetry of H. For any operator 0{v) 
involving only gravitational and fields, -^{0{X,T)) can be calculated using H{Np)\^) = to 
obtain 



d_ 
df 



{0{X,T)) = {^\ 



Tl 



d{X),H[N;i 



(4.8) 



where the 0{h) correction comes from the quantum fluctuation of the clock fielqj [48j. To carry 
on, we first recall that our specific matter coordinates and frames satisfy (3.30). The correspon- 
dence (3.34) implies that under a change of the spatial coordinates and frames, the emergent fields 
{E'}{X,T),A({X,T),(t){X,T),P{X,T)] transform classically by passive local SU{2) and diff m 
transformations [H], up to errors of 0{h) + 0{d). Therefore, by applying the classical transfor- 
mations to the result of our calculation, we can generalize it to arbitrary spatial coordinates and 
frames, up to errors of 0{h) + 0{d). 

In the matter coordinates and frames satisfying (3.30), the applications of (4.8) lead to 



d 

dT 



d_ 
df 

m 



EnT)dsa 
Tl JSx,Ax 

^ Jiex,Ax)i,Hg^^ (Npi 



d 

dT 



Tl 



d 

dT 

{^\\ 



Vexpif Ai{T)irj)de% 

J ex r\x 



h{ex, 



ex, AX 
Ax)j 



N } 
p* 



liTA"^) + e{^} + o{n) 



h{exAx)j,'Hg,^ (i^i HtiI^) + e{^} + 0{h) 



(4.9) 



_d_ 
dT 



{piX,T,)) 



Tl 



d 

dT 



I P{X',T)dX' 

Tl Jcx 



N i 



tlrM + e 



0{h) 



{^^l [p(X),7i,,^ (z.f 1)1 IiTi\^) + e{^} + 0{h) 



where we again replace Hg^^{Np*") with 'Hg^^{Afn) using the spatial coordinate condition (3.21). 

^In general, there is also a correction of 0{d) in (4.8) due to the discretization error of inhomogeneous 
VijjiX, Tl). In this paper, we have a homogeneous i/^ field at Ti, so the correction does not appear here. 
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We then apply (3. 28), (3. 19) and (3.31) to obtain: 



Ti 



J 



i{J{T,)MTi),h'^'\T,),p{n)m 



J ) 



i{j{n)),{m)),Ch(>^\T,)),{p{Ti))) 



d 

dT 



X,AX\K 
h 



'[{J{Ti),h{Ti),h<^^^iTi),p{Ti)m + e{^} + 0{h) 
= (^J'^'')f mTi)), {h{Ti)), {h^'HT,)), {p{T,))) + + 0{h) 



{p{x,Ti)) = p{x),ng^^o Ko-^(ri)) HtJ*) + + o(^) 

= (*|$^(J(Ti), /i(Ti), h(^\T^),p{T,)m + e|^} + ©(fi) 

= $^((J(ri)), (/i(ri)), (/iWm)), (p(Ti))) + e^^y + o(/i) 



(4.10) 

The first terms in (4.10) are the 0(hP) contributions ignoring the back-reaction of {ip} sec- 
tor. They are given by the functions ($7 )/, (^'j,^ )/' iK )l and which are func- 
tions of the emergent fields through (3.34). Recall that in our spatial coordinates AXj„^j G 
{(ltd, 0, 0), (0, ltd, 0), (0, 0, id)} with ^ {^n^}' and define (5^^ to have the only none-zero com- 
ponents: 

s{±dm ^ ^ ^(o,±d,o) ^ ^1 ^ ^(o,o,±d) ^ 

For our symmetric case (4.6) we have: 
($J'^^),((J(ri)), (/i(ri)), (/iW(Ti)), (:p(ri))) = 5f^$j(d2£;i,exp(d^ir,),exp(U0i),cZ3Pi) 

{^ia^'mhTi)), {h{Ti)), {h^'Hn)), {p{T,)), (a(ri))) = 

(^f ^'')f ((-^m))' (^^^H^i)), {p{Ti))) = St''{r')fMd''E^,exp{dA,Ti),eMi>^h),d^Pl) 



^^i{J{Ti)), {h{Ti)), (/iW(ri)), {piTi))) = ^p{d''Ei,eMdAiTi),eMiMi),d^Pi) 



(4.11) 



where the new set of functions $p} have no labels X, AX or /. The term vanishes 

because of the diagonal condition on Ef in (4.5). Refering to (4.10) and (4.9), we see that (4.11) 
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implies the preservation (4.6) along the clock time. Thus, at 0{hP) and ignoring the back reactions 
from m', the dynamics of the gravitational and ^ field is isotropic, homogeneous and spatially flat: 



Ef{X, T) = 5'iE{T) + £{^1 + 0{h) ■ AiiX, T) ee 5iA{T) + e|^} + 0{h) 
P{X, T) = P{T) + e{^} + 0{h) ■ <t>{X, T) = cf>{T) + e^^} + 0{h) 



(4.12) 



with the initial condition (4.6): E{Ti) = Ei, A{Ti) = Ai, P{Ti) = Pi and 0(ri) = 

Recall from (4.5), the emergent fields satisfy the Gauss, momentum and Hamiltonian constraints 
up to the corrections. In our specific case, the Gauss and momentum constraints are trivially 
satisfied while the Hamiltonian constraint gives nontrivial implications for the dynamics. Inserting 
(4.12) into (4.4) with n = 1, we have: 



<^i{<fE{T), eMdA{T)rj),eMi^<l>{T)),d^P{T),M) + e^^y + 0{h) 

= Hg,<t>i^)\siE{T),S(A{T),(j>(T),P{T),N=M + ^(^^) + + ^(^) 

= 



(4.13) 



Here we see that the 0{hP) effect ive Hamiltonian constraint for gravitational and ^ field is given 
by $1, which contains corrections of O(ci^) to the classical -ffp,,/,. 

To describe this symmetric sector in a symmetrically reduced form, we introduce the sym- 
metrically reduced variables - {p,hi = exp(crj)} {i = 1,2,3) in place of {cP E , exp^dArj)} , and 
{p^,h(j, = exp(zAo0)} in place of {(i^P,exp(iAo<^)}. Their non-zero Poisson brackets are defined 
by: ^ 

{p,c}=-Kj- {p^,0} = l (4.14) 

where the factor | accounts for the degeneracy on the three independent spatial directions. We 
now define Hiqc{N) {N G M) to be a function of {p, hi,p^, h^} that satisfies: 

^/.c(iV)|^^,,^,^^^^^p(,^^^)_^^^,3p,^^=exp(ao.^) = *i(ci'£?,exp(dAr^),exp(ao0),d^P, JV) (4.15) 
Its explicit form is given by: 
H,„(N) = if« (JV) 



hr\h,,{H^{l),p^}hJ'{hJ,{H^{l),p^}h,{h,\p^ 



-1 



-1 



i,j,k=l 



+N ■ 2Kp-3/2p2 



where 



Next, we define: 



i,j,k=l 



0" 



2iX( 



-[h^-hl\ - t;^o = {0O,i/,,e(l)} 



(4.16) 
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and find that the equations of motion (4.12) can be expressed as 
d^E{T) = {p,Hig,{N)} 
+e{^} + 0{h) 



dT 



N=vJ; p=d^E(Ti),hi=eMdA{Ti)Ti},p^=d3p{n),h^(Ti)=eMi^o<f>(Ti)) 



_d_ 



eMdA{T)n) = {hi,Higc{N)} 
+e{^} + 0{h) 



N=vJ; p=d2E{Ti),h.i=eMdA{n)Ti),p^=d3P{Ti),h^{Ti)=exp{i\o(i>{Ti)) 



d_ 



d'P{T) = {p^,Hi,,{N)} 

I 

+e{^} + 0{h) 



N=vJ; p=d^E{n),hi=eM<iA{Ti)Ti},p^=d3p{n),h^(Ti)=exp{iXo<i>{Ti)) 



(4.17) 

Equations (4.17) state that the 0{h^) evolution of the emergent gravitational and (f) fields in the 
state 1^') is governed by the effective Hamiltonian constraint Higf.{N), ignoring matter back reac- 
tions from {'(/)} sector. To this approximation, the symmetrically reduced model captures the clock 
time dynamics of the emergent fields through (p((^°), c(<^°),p^(<^°)) = {d^E{T),dA{T),d^P{T)), 
and it is straightforward to check that we have a conserved observable d^P{T) = const = d^Pi. 

Next, we investigate the evolution of the spatial scale. Recall from (2.16) that the volume 
operator with a dynamical region R is obtained by summing over the operators with Vn G R- 

Using (3.15) and (3.17), we can localize \/\Q{Xn, Therefore, the spatial volume 

observable corresponding to the coordinate region Cl C P at a clock time T is given by: 

vin,T)= J2 V\QiXm,T)\ 

To investigate the evolving spatial scale in the cosmology, we now keep track of the volume of the 
spatial region coordinatized by cx^ at various clock times: 



d_ 
df 



Ti 



d_ 
df 
_d_ 
df 



e^JKeabcE'}E^jE1^{X',T)d^X' + 0{h) 



d^E^/'\T) + 0{h) 



N=vJ; p=d'^E{Ti),hi=eMdA{Ti)n),p^=dSp{Ti),h4Ti)=eMi>^0't>{Ti)) 



+e{^} + o{n) 

-o ^^7f^d^E{Ti) [2sin(o!A(ri)cos(dA(ri) [l - 2(1 + 7^) sin2((iA(ri)]] 

+e{^} + 0{h) 



(4.18) 
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Also, the constraint equation (4.13) leads to: 



(4.19) 



Since the energy density of the (f) field is given by p^{T) = Pc is equal to the conserved 

quantity {§jK^-i^)-^/''{d^P{T))-^ = {§^KS^)-y^{d^Pi)-\ 

To obtain our modified first Priedmann equation, we solve the constraint equation (4.19) and 
find: 

2 

sin2(dA(T,)) = ' +"7^'^ + + ' " + (^) ' ^'-'^^ 

Substituting (4.20) into (4.18) we obtain the first modified Priedmann equation for the Hubble 
constant H{T): 

(4.21) 

When the universe is in the classical region, we have ^ P</>(^i) and x(Ti) ^ 1, and (4.21) 
approaches the classical first Priedmann equation: 

H\n) = |p^(Ti) + e{^} + 0{h) (4.22) 

When becomes comparable to p^, our model gives significant modifications to the classical 
FRW cosmology. It is clear from (4.21) that H{Ti) vanishes when p^(Ti) = p^ and that the initial 
singularity is replaced by a bouncing behavior, up to the corrections of e^^y + 0{h). 



5 Comparisons with Loop Quantum Cosmology at 0{ff) 

We have now identified a symmetric sector of the model described in the first half of the paper. 
This sector is represented by a state l^*) which gives an approximately homogeneous, isotropic 
and spatially flat evolution of the emergent gravitational and (p fields. Further, the evolution is 
effectively governed by the symmetrically reduced classical Hamiltonian constraint Hiqg{N). Also, 
the evolution of the fields agrees with PRW cosmology in large scales, while it deviates from PRW 
cosmology when the cp field energy density is close to the critical value p^.. Finally, the deviation 
results in the resolution of the initial singularity. In the following, we will compare Hiqg[N) with 
the effective Hamiltonian constraints of 0(hP) in the existing models of loop quantum cosmology. 

The homogeneous, isotropic and spatially fiat sector of classical general relativity can be de- 
scribed by a symmetrically reduced theory, namely the flat FRW cosmology. In our setting, the 
flat FRW cosmology describes an arbitrary representative cell C of the comoving space with four 
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(off-shell) phase space variables {p,c,cf),p^}. Since the space is flat, we can choose an arbitrary 
Euclidean spatial coordinate system for C, in which C has a coordinate volume V. In this spatial 
coordinate system, the reduced phase space variables are related to the Ashtekar's variables by 
[271 [281 [29]: 

V'/^Ai{x) = c-ooi; V'"^E'j{^)=p-e]-, VP(x) = p^; (/>(x) = (5.1) 

The triads and cotriads are non-dynamical, and they will be gauge-fixed to identity matrices 
from now on. The two conjugate pairs satisfy the Poisson brackets: 

{c,p} = ^K7; {</>,P0} = l 
For our case, the Hamiltonian constraint HpRwi^) for flat FRW cosmology is given by [27 1 [28 t [29] : 

Hfrw{N) = -N ■ -^c'p^/^ + N • 2Kp-3/2p2 (5,2) 

where iV G M is the symmetrically reduced lapse function. 

In the spirit of loop quantum gravity, loop quantum cosmology uses {p, /if = exp(/icrj)} as 
elementary gravitational variables instead of {p,c}. The holonomy is given by the parallel 
transportation by the connection field A^^ along a certain path in the direction of e". In the Eu- 
clidean spatial coordinate system used in (5.1), the coordinate length of the path is set to be fiV^^^. 
Each model in homogeneous, isotropic and spatially flat loop quantum cosmology has an effec- 
tive Hamiltonian constraint of 0{hP) that approximates Hfrw{^) in terms of the loop variables 
{p, /if}. Clearly, there are ambiguities in such approximations, and the different approximation 
schemes based on different physical considerations result to the different models in loop quantum 
cosmology [M ISSl [26l [IB [2S1 [291 [30j • 

The first main ambiguity lies in the ordering of the two required procedures in obtaining the ef- 
fective Hamiltonian constraints of loop quantum cosmology from the full form ff^ ,^(iV): symmetry 
reduction of the phase space and the introduction of the loop variables. When the symmetry reduc- 
tion is applied first [271 [281 [29j , Hg^^{N) is first simplified to the gravitational term in Hpr\y{N), 
before the replacement of {p,c} by {p, /if}. When the loop variables are introduced first [30] . 
Hg^(f,{N) is first approximated by the local holonomy and flux variables, before imposing the sym- 
metry on the local loop variables and reducing them into the set {p, /if}. The first procedure leads 
to a simpler effective constraint, while the second leads to an effective constraint closer to the full 
form. We will denote the two schemes as / and //. 

The second main ambiguity lies in the choice of C and ^. In FRW cosmology, the choice of 
C does not effect the physical predictions for {^^, P}. This is because of the consistent 

scaling of {p,c^cf),p^} and HpRwi^) when C is changed. Given the scaling of {p,c,(f),p^} with 
V according to (5.1), HpRwi^) scales linearly with V as supposed. In loop quantum cosmology, 
this independence from the choice of C may be disrupted due to the replacement of of c by /if . 
Particularly, the effective hamiltonian constraints of loop quantum cosmology can scale nonlinearly 
with V when C is changed. Therefore, the choice of C does affect the physics in some models loop 
quantum cosmology. 

Given a choice of C, the value of fi must also be specified. Clearly, the value of ^V^/^ must be 
small enough that the effective Hamiltonian constraint can give a good approximation of HpRW- 
In the early stage [25j of loop quantum cosmology /i is set to be a small fixed constant fi = fiQ. 
However, two main objections were raised against this choice [M]- First, with /x = fiQ, the 0{hP) 
effective Hamiltonian constraint scales nonlinearly with V under a change of C. This is a problem 
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for loop quantum cosmology since the physics of {A^^, Ej, (p, P} depends on the choice of C, while 
there is no natural way of making such a choice. Second, with /i = ^uq, the critical (j) held energy 
density pc at the big bounce is inversely proportional to the conserved quantity P. Thus there is a 
danger of predicting the big bounce at a low matter density if P is high enough. This is a problem 
since a high P value is preferred for semi-classical limits. Because of these objections, a physically 
motivated choice with ii = Ji = Dp^^/"^ has been adopted [211 |26l |27], where D is a real fixed 
parameter. The appearance of p in restores the independence of C for loop quantum cosmology 
with ^ = p.. Moreover, it also leads to a fixed critical density set at the Plank scale. We will denote 
the old and new choices in p and C as (C,/xo) and (C,/i) schemes. 

The combinations of these schemes give four distinct versions of the 0{hP) effective hamiltonian 
constraint H^^^. We will use the superscripts {C,fj,o,I), {C,fio,II), {C,fl,I) and {C,p,,II) to 
indicate the schemes applied. First, we define the H^^^ with a general choice of /u as 



i,j,k=l 



{K-h%)hl'^{hlpi] 



The two effective Hamiltonian constraints using fi = fiQ are: 



(5.3) 



LQC 



47V(l + 72) 



4^,7 



(5.4) 



The other two using fj, = fj, = Dp 



-1/2. 



Hf£'\N) = -r^H^L^^) + N2kp''/'pI 



4 4iV(l + 72 



9 



-iV2Kp-3/2p2 



t,j,k=l 



One can easily check that 



rE,il 



(5.5) 



HFRwiN)= lim HffJ>''\N)= hm HffJ'''\N) = lim Hf'^f^^'\N) 



(5.6) 



lim H'^^'P'^^Kn) 



26 



Note that (5.3) and (5.4) have almost the same form except for the extra 4/9 factor in H^^^^'^^\ 
and their quahtative difference hes in that p, depends on p. 

The simplest (C, hq, I) scheme was studied in the early stages [25] of the program. The improved 
scheme {C,fi,I) has been extensively studied both analytically and numerically [261 123 |28l [29] . 
The {C,fL,II) scheme that is closer to the full form of loop quantum gravity has recently been 
investigated Lastly, the effective model we derived from the 0{hP) contributions of the semi- 
classical limit of the full theory corresponds to the (C, /io, //) scheme. Denoting cx^ as the region 
of the emergent space coordinatized by cx^, we compare (4.15) and (4.16) to (5.4) and find: 

^£r'''^(^)lc=c,„,,o=i = ^"^^(^) (5-7) 

While the different schemes are physically distinct, they share important qualitative features, 
which represent the most robust part of loop quantum cosmology. The first three schemes consis- 
tently resolve initial singularity by a bouncing of the scale factor, as shown in [2^1 [271 [2EI [221 [30] 
to various approximations. In the approximation of 0{hP), we have seen that the semi-classical 
limit given by the state |^') in our mode matches with the {C,ij,q,II) scheme, which also predicts 
the same bouncing behavior. From (4.20) and (4.21) it is clear that when the density is close 
to the critical density p^, the contributions from sm{dA(T)) that are nonlinear in dA(T) become 
important, and drive the evolution away from the FRW cosmology. At 0{hP) level, the big bounces 
in the four schemes are caused by these non-linear holonomy corrections. From (5.6) we see that the 
holonomy corrections appears only with non-zero values of fiV^^^, which via (5.7) corresponds to 
the non-zero value for d in our model. Since the non-vanishing d in our model reflects the intrinsic 
discreteness of space, our model supports the idea that the quantum geometry leads to the big 
bounce. 

It should be emphasized that we obtained the {C,fio,II) scheme from the spatial quantum 
geometry of loop quantum gravity. Through our model, the first-principle derivation not only leads 
to the specific scheme, but also determines the values of C and fiQ with fundamental meanings. To 
see this, let us look at the objections to the /iq schemes in light of our model. Recall that the first 
objection involes the sensitivity of the fiQ schemes to the arbitrary representative cell C in the space. 
From (5.7), it is clear that such sensitivity corresponds to our model's sensitivity to cx„- However, 
the objection does not apply to our model, where cx„ is not arbitrary, but is an elementary cell in 
the space that is dual to a single physical node. The elementary cell represents one of the smallest 
regions in the emergent space with nonzero volume. Conversely, our model states that the cell C 
in the /Uq scheme models corresponds to one of the elementary cells of the emergent space in the 
full theory, and that /io = 1 so the holonomy variables simply run over the sides of this elementary 
cell. 

The second objection involves the possibility of the big bounce happening at a low matter 
density, or at a large volume of space. In the context of our model, one can explicitly check 
whether this really occurs for the state l^*) describing our universe. Through the model, we now 
make a trial estimation of the scale of the big bounce using known cosmological factors. Recall 
from (4.19), (4.20) and (4.21) that the (p field energy density is p^{T) = -^sjj^ and the critical 
energy density is p^ = (^K'^7^)~^/^((i'^Pi)~^. At the moment of the big bounce we have: 

On the other hand, at current cosmological time Ti we have: 

, , 2K(d3Pi)2 ^ ^ 
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In our setting with negligible energy density of the {V'} sector, p^{Ti) is approximately equal to 
the dark energy density observed today, which is p^(Ti) ~ 10~^ J/m^. Here we will use the value of 
Immirzi parameter 7 ~ ln(2)/7r, given by black hole entropy consideration in loop quantum gravity 
|47j . Also, since no dispersion anomaly has been detected even in the cosmic-ray protons with 
wave lengths of 10~ m, we expect £E{Ti) < (10-26)2 

so the space appears smooth for those 
protons. Applying these values to (5.8) and (5.9), one finds the volume of cx^ at the big bounce 
to be d?E'^/'^{Tc) « lO^^^^m^ compared with the Planck volume Vp lO^^^^m^. Therefore, to the 
0{hP) approximation, the big bounce does not happen at the large volume region as one might 
worry. On the contrary, the higher order contributions in h for this case are obviously important 
near the big bounce given by the 0{hP) contribution. 

The source of these higher order contributions contains not only quantum fluctuations of the 
scale factor, but also inhomogeneous and anisotropic quantum fluctuations. Clearly, we may cal- 
culate these corrections only if we explicitly construct the coherent state | ^) . This is an important 
task for the development of the model, and hopefully it will yield further insights for loop quantum 
cosmology. 



6 Conclusion 

This paper starts from the kinematical Hilbert space of loop quantum gravity, which describes 
the matter fields living in the dynamical quantum geometry of space. Using the model with a 
modified Hamiltonian constraint operator, we see that the dynamics of such a system reproduces 
FRW cosmology in the large scale limit. Further, the 0{h^) corrections of the model for FRW 
cosmology conform with loop quantum cosmology in a specific scheme. Such a result is valuable, 
since it attributes the predictions of loop quantum cosmology to the fundamental principles in loop 
quantum gravity. 

The result serves as a starting point to many possible future projects. First, one may explicitly 
construct the coherent states in the model to evaluate the emergent cosmology beyond 0{hP), to get 
the quantum fluctuation corrections in the emergent cosmology. Second, one may try to derive more 
of the implications of loop quantum cosmology by applying the model to more realistic cosmological 
settings. Third, one may try to improve the model by incorporating the graph-topology changing 
feature in the Hamiltonian constraint operator, in the hope of deriving loop quantum cosmological 
models with /J, = fi. 
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